A0-A119  10S  SCORStA  INST  Of  TCCH  ATLANTA  SCHOOL  OF  PHYSICS 


UNCLASSIFIED 


CORRELATION  and  collective  noocs  in  narrow  SANO  NATERIALS.(U) 
HAT  SI  MW  RIBARSKY  AFOSR-AO-0023 

AFOSR-TR-BB-OTll  I 


F/8  10/14 


N 


AD  A119188 


iwwiaaiBi 


JtCURITY  tl  ASM F|C  ATION  O*  i  Ml*  HKUt  (*»n*n  umt*  e.m«r«oj 


REPORT  DOCUMENTATION  PAGE  beforeVompletwc^orm 


t.  REPORT  NUMBER  .  |2.  GOVT  ACO'.ittOM  NO.  ».  RECIPIENT'S  CATALOG  NUMBER 

AFOSR-TR-  82-071  1  | ftD-A//9/2$ 


4.  TITLE  (and  SuAilll*) 

"Correlation  and  Collective  Modes  in  Narrow  1  / 
Band  Materials". 


7.  author^*; 


Martin  W.  Ribarsky 


».  type  of  report  a  pcrioo  covcreo 
Interim  1980 

Oct  1,  1980  to  Sept. 30, 1981 


s.  performing  org.  report  number 


1.  CONTRACT  OR  GRANT  NUMBERS) 


AF0SR-80-0023 


S.  PERFORMING  ORGANI2ATION  NAME  ANO  AOORESS 

Georgia  Institute  of  Technology 
School  of  Physics,  Atlanta,  Ga.  30332 


11.  CONTROLLING  OFFICE  NAME  AND  AOORESS 

Air  Force  Office  of  Scientific  Research 
Bolling  AFB,  Bldg.  410 

Washington,  D.C.  20332 _ 

4.  MONITORING  AGENCY  NAME  ft  ADORESSI'I# S«MmI  Horn  "ConUmtUoi  Of  tic*) 


16.  DISTRIBUTION  STATEMENT  ( ©/  thlm  Report) 

Approved  for  public  release; 
Distribution  unlimited. 


10,  PROGRAM  ELEMENT,  PROJECT.  TASK 
AREA  ft  PORK  UNIT  NUMBERS 

fcfl  02  P 
A30|  /A* 


«.  REPORT  OATt 

May  25.  1982 


IS.  NUMBER  OF  PAGES 


IS.  SECURITY  CLASS,  (ml  Nil*  report) 


iCUSSiflED 


DECLASSIFICATION/  DOWNGRADING 
SChCOULE 


r VQ 


If.  KEY  RORDS  (Contlnum  on  rmvmttm  aid*  It  nmcmmmmry  mod  Idmntlfy  by  block  ntmkmr) 


Transition  metal  oxides,  correlation,  excitons,  interband  excitations,  energy 
loss  spectra,  superconductivity,  band  structure,  cuprous  oxide. 


20.  ABSTRACT  (Contlnum  on  rmvmrmm  mi  dm  II  nmcmmmary  mnd  Identity  by  block  rutmbmt) 

^Correlation  and  collective  modes  have  been  studied  for  systems  with  quite 
localized  valence  or  conduction  bands.  In  particular  this  research  has  been 
concerned  with  localized  electron-hole  states  and  how  they  contribute  with 
other  excitations  to  the  dynamical  response  of  the  system.  The  most  important 
aspect  studied  has  been  the  effect  of  excitons  or  exciton-like  states  on 
superconducting  properties.  Initially  the  system  studied  has  been  Cu^O  for 
which  a  simplified  tight^binding  model  was  used.  The  results  showed  that  an 


KWiil 


effective  attractive  interaction  between  electrons  due  to  the  local  field 
'  produced  by  interacting  electron-hole  pairs  was  possible.  Under  the 
appropriate  conditions  this  interaction  could  lead  to  a  superconducting 
state.  The  contribution  of  the  localized  excitation  structure  and 
plasmon  modes  to  the  electronic  loss  spectra  was  also  analyzed.  Here  the 
effect  of  the  primary  electron  energy  on  the  plasmon  peak  position  was 
shown  to  be  significant. 


Accession  For^ 

Intis  •  gra&i 

DTIC  TAB 

Unannounced 

Justification — 


By - - - - - 

Distribution/ _ 

Availability  Codes 
lAvail  and/or 
>  Dist  !  Special 


mm j 


SECURITY  CLASSIFICATION  OF  THIS  PAGEr**«n  Omtm  Rnfr*4) 


Table  of  Contents 

Page 


1.  Research  Objectives . . . .  1 

II.  Localized  Effects  in  Superconducting  Systems.... ...... .................  3 

A.  Alternative  Coupling  Mechanism... . . . .  3 

B.  Excitonic  Effects.... . . . . . 5 

C.  A  Tight-Binding  Approach  to  Calculations  of  Excitation  Structure...  8 

III.  Development  of  Calculation  Methods... . . . 13 

A.  Molecular  Cluster  Methods . . . 13 

B.  Band  Structure  Methods . . . . . . .  14 

IV.  Electron  Energy  Loss  Structure  of  Clean  and  Oxidized  Titanium..........  15 

V.  References............. . . 18 

VI.  Publications  and  Presentations . 20 

VII.  Personnel .  21 

VIII.  Interactions... . . . . . 22 


localized  collective  modes  to  understand  the  interaction  with  optical  or  elec¬ 
tronic  probes  or  large  scale  correlated  states  such  as  superconductivity*  In 
certain  cases  one  must  consider  these  effects  to  understand  even  ground  states 
properties,  such  as  for  narrow  band  insulators*  In  particular  this  research  has 
been  concerned  with  localized  electron  -  hole  states  and  excltonic  modes  and  how 


they  contribute  with  other  excitations  to  the  dynamical  response  of  the  system* 
During  the  preceding  year  we  have  completed  studies  of  interband  excitations 
and  surface  and  bulk  plasmon  modes  and  their  contributions  to  the  electron 


energy  loss  spectra  of  TiO  and  FeO.  These  studies  were  important  for  understanding 
the  structure  of  the  loss  spectra,  which  can  be  a  very  precise  probe  of  surface 
electronic  structure  and  the  effects  of  adsorption. 


We  have  also  begun  studies  of  mechanisms  for  superconductivity  which  are 
alternatives  to  the  phonon-coupling  mechanism*  Among  these  new  mechanisms  are 
coupling  via  plasmons  or  excitons*  The  excltonic  mechanism  is  especially  promising 
for  narrow  band  systems,  and  we  have  started  investigations  of  Cu^O  to  ascertain 
,  whether  excltonic  effects  can  cause  an  attractive,  non-local  dynamical  response* 
This  work  may  lead  to  an  explanation  of  some  unusual  experimental  results  for 
this  and  other  systems* 


In  conjunction  with  these  investigation  of  particular  effects  for  specific 
systems,  we  have  continued  a  program  of  developing  computer  codas  useful  in  the 
study  of  localized  systems*  These  include  band  structute  codes  and  cluster  codes 
for  studying  well-localized  structures  at  the  surface  or  in  the  bulk*  With  the 
band  structure  codes,  ve  have  expanded  the  number  of  lattice  symmetries  available 


(1) 


to  Include  the  rocksalt  structure  appropriate  for  transition  metal  monoxides • 
Vlth  the  cluster  codes*  we  have  added  new  features  to  Improve  convergence  and 
to  simulate  better  the  bulk  environment  for  clusters  representing  defects  in 
the  solid. 

In  the  following  sections,  we  will  discuss  how  we  proceeded  toward  our 
objectives.  Ve  will  also  discuss  some  of  the  future  work  we  will  do  in  this 
area. 


(2) 


II*  Localized  Effects  In  Superconducting  Systems 


A*  Alternative  Coupling  Mechanisms 

From  time  to  time  some  researches  have  proposed  and  others  have  looked  for 
other  mechanisms  besides  the  basic  phonon  mechanism  for  obtaining  the  effective 
attractive  interaction  between  electrons*  The  search  was  motivated  in  part  by 
the  enormous  benefit  which  would  accrue  from  finding  a  superconducting  state  which 
persisted  even  up  to  liquid  nitrogen  temperatures  (around  70°K) .  The  highest 

critical  temperature  yet  found  is  in  the  range  of  25-30°K* 

1  * 

Bardeen,  Cooper  and  Schrieffer  showed  that  superconductivity  arises  from 

an  effective  attractive  interaction  between  quasi-particles  near  the  Fermi  surface 
of  a  degenerate  electron  system*  In  their  work  and  in  most  subsequent  work,  the 
effective  interaction  is  assumed  to  be  a  constant  over  frequency  up  to  a  cut-off 
frequency  beyond  which  it  is  zero*  However,  the  effective  interaction  need  not 
be  constant  and  it  need  not  be  due  to  phonon  coupling  alone*  Some  worker  have 
considered  additional  mechanisms  such  as  plasmons  or  excitons  and  have  also  con¬ 
sidered  a  dynamic,  momentum-dependent  interaction*  It  has  been  speculated  that 
2  3  4 

plasmon  or  exciton  *  coupling  can  greatly  increase  the  critical  temperature  Tc* 

Different  structures  such  as  layered  metal-oxide-semiconductor  systems  or  thin 

films  have  also  been  considered  as  candidates  for  enhanced  Tc* 

An  electron-hole  plasmon  mechanism  was  shown  by  Pashitkli^  to  give  an 

attractive  coupling  for  model  semiconductor  or  semi-metal  systems  with  appreciably 

2 

different  conduction  and  hole  electron  masses*  Later  Tanaka  showed  that  the 
electron  plasmons  alone  would  support  the  attractive  coupling  for  semiconductors 
In  both  two-  and  three-  dimensional  systems*  However,  Kelly  and  Hanke^  did  a 
careful  calculation  for  Sl(lll)-  SiO^  in  a  MOS  system  where  they  considered  a  coupled 
electron-phonon  and  electron-plasmon  mechanism*  They  solved  the  veak-coupling  super¬ 
conductivity  equations  for  a  number  of  carrier  densities  but  found  only  very  slight 

•  \ 

enhancements  of  T  due  to  plasmons  at  low  carrier  densities*  No  enhancement  of  T  i 

c  c 

which  could  be  attributed  to  plasmon  coupling  has  been  seen  experimentally.  j 


Other  coupling  mechanisms  have  been  considered  which  are  related  to  the 
plasmun  or  excitonic  mechanisms*  An  acoustic  plasmon  mode  which  can  be  due  to 
Interacting  electrons  and  holes  as  in  the  Fashltskii  model^  or  to  the  interaction 

between  heavy  d  and  light  s  electrons  has  been  suggested*  Calculations  for  tran¬ 
sition  metals  using  the  latter  interaction  show  that  the  enhancement  of  T  is 

c 

minor Another  model  involves  an  electron  heavy  hole  mechanism  where  the  electron 
pairing  arises  from  the  exchange  of  virtual  bosons  associated  with  vibrations  of 

g 

the  hole  lattice*  Here,  again,  more  complete  calculations  show  that  the  enhancement 
o 

of  T  would  be  small  * 
c 

3 

The  excitonic  mechanism  of  superconductivity  has  been  considered  before  ,  but 

4 

the  most  recent  model  is  due  to  Ailender  Bray  and  Bardeen  (ABB).  They  chose  to 
consider  excitons  in  a  semiconductor  where  they  are  well-defined  excitations.  A 
metal  in  Intimate  contact  with  the  semiconductor  provided  free  carries  which  inter¬ 
acted  by  way  of  the  excitons.  The  authors  chose  a  narrow-gap  semiconductor  to 
treat  within  their  model  since  the  valence  electrons  can,  in  many  cases,  be  con¬ 
sidered  as  if  in  a  free  electron  metal  with  the  gaps  introduced  by  perturbation 
theory.  The  model  calculations  were  made  to  second  order  in  the  electron-exclton 
coupling  and  revealed  an  attractive  local  polarization  field  due  to  the  exciton 
coupling.  Under  certain  very  favorable  conditions  the  model  predicted  a  value  of 

T  at  room  temperature  or  above*  However, a  later  calculation  by  Bardeen^  in- 
c 

dlcated  a  much  smaller  enhancement  of  T  was  likely. 

c 

Inkson  and  Anderson  later  came  out  with  a  rebuttal  to  the  ABB  arguments.  ^ 

They  showed  that  the  use  of  second-order  perturbation  theory  could  lead  to  double 

counting  of  the  exciton  response  and  that  one  should  use  the  complete  inverse 

dielectric  function  e~^  for  the  semiconductor.  Their  calculations  using  a  long- 

wavelength  model  dielectric  function  for  silicon  and  indium  antimonide  showed  no 

enchancement  due  to  excitons,  as  was  also  the  case  with  later,  more  detailed  cal- 

12 

culations  for  Ge  using  pseudopotential  energy  band  structure  and  wavefunctlona . 

13 

ABB  countered  that  the  Inkson  and  Anderson  approach  was  valid  for  the  systems 
they  considered,  but  that  the  ABB  model  was  for  different  systems*  ABB  said  that 
excitonic  effects  were  more  likely  in  systems  where  excitonic  excitations  were 


strong  and  local  field  effects,  which  Inkson  and  Anderson  had  neglected,  were 
Important* 

There  is  some  experimental  evidence  which  may  support  the  ABB  contention* 

Researchers  have  postulated  a  similar  excitonic-superconductivity  mechanism  is 

responsible  for  recent  observations  of  near-perfect  diamagnetism  in  temperature- 

o  14 

quenched  CuCt  and  pressure-quenched  CdS  at  temperatures  up  to  300  K  •  Measure- 

15 

ments  by  Lockhart  etal  have  established  the  existence  of  an  anomalously  weak 
electric  field  outside  the  surface  of  copper*  Below  4.5°K  a  phase  transition 
occurs  and  the  electric  field  drops  from  the  expected  value  to  near  zero.  A 
thermal  oxide  some  20  X  thick  covers  the  Cu  surface,  and  the  perfect  shielding 
most  likely  occurs  in  that  region.  The  most  likely  explanation  of  the  shielding 
is  that  the  surface  electrons  couple  and  form  a  Bose  condensate. ^  It  is  also 
possible  that  this  phase  is  superconducting.  Other  experiments  indicate  that  an 
exciton  gas  formed  in  Cu^O  by  laser  excitation  obeys  Bose-Einstein  statistics.^2 
Since  Cu^O  is  a  semiconducting  material  with  strong  excitonic  effects  of  the  type 
considered  by  ABB,  the  anomalous  behavior  may  be  due  to  their  mechanism. 

B*  Excitonic  Effects 

The  usual  effective  interaction  constant  for  superconductivity  is  written 

N(0)V-  y-X  (1) 

where  N(0)  is  the  density  of  states  at  the  Fermi  level  of  one  spin,  y  is  the 
density  of  states  times  an  averaged  screened  Coulomb  interaction  and  X  is  the 
attractive  coupling  constant.  The  effective  interaction  may  be  written 

V(q,«)»  4lle2/q2e(q,w) 

for  vavevector  q  and  frequency  u  where  c  is  the  total  dielectric  function  con¬ 
taining  both  electronic  and  Ionic  polarizabilities.  Now  stability  requires 
that  c(q,o)  >0  so  by  Eq.(l)A  <  |i  at  least  at  u  «  0.  However  it  is  known  that  X 
can  be  auch  larger  than  p  in  many  cases.  Cohen  and  Anderson  have  resolved  this 
contradiction  by  pointing  out  that  the  stability  requirement  applies  only  to  the 

ionic  interaction  and  that  the  polarization' fields  seen  by  the  ions  and  the 

(5) 


electrons  are  different,  the  difference  being  that  the  electrons  see  the  ionic 

18 

self-polarization.  When  ionic  local  field  effects  are  included  ,  one  gets 
2  2 

X«  / c^h  where  is  the  ionic  plasma  frequency  and  is  a  phonon  frequency. 


Now  X  can  be  larger  than  p,  and  since  the  McMillan  result  for  T  is 

c 


19 


T  -  expf  _ I 

c  expl  A-  ji(1+.  62A)  J 


(2) 


one  can  have  higher  transition  temperatures. 

4 

Allender,  Bray  and  Bardeen  postulated  that  a  similar  local  field  arose 

2  2 

for  the  exciton  mechanism  and  that  for  a  semiconductor,  X  «  m  /w  where  w  is 

ex  p  g  p 

the  plasmon  energy  and  u)  is  the  gap  energy.  They  derived  X  for  a  narrow  gap 

g  ex 

semiconductor  and  showed  that  in  favorable  cases  with  their  model,  one  could  get 


X  >  p  and  substantial  increases  in  T  . 
ex  —  c 

As  was  stated  in  Sec.  II.  A,  the  main  objection  of  Inkson  and  Anderson 


11 


to  the  ABB  model  was  that  one  should  use  the  complete  inverse  dielectric  function 
e  *  rather  than  perturbation  theory  to  calculate  X^.  They  also  said  that  the 
exciton  and  other  excitations  are  poles  of  the  dielectric  function,  but  are  zeroes 
of  the  interaction  because  that  involves  1/e.  Thus  th$y  argued  that  excitons 
should  not  be  important  to  the  effective  interaction  and  backed  this  contention 
up  with  model  calculations  of  e  *  for  Si  and  In  Sb.^  One  can  show,  however,  that 
this  last  argument  of  Inkson  and  Anderson  is  not  true  in  all  cases.  To  see  this 
we  consider  a  system  with  an  Insulating  ground  state  in  the  Wannier  representation. 
The  local  Wannier  functions  and  energies  are  related  to  their  counterparts  in  the 
Bloch-representatlon  by 


„„<?-$  -  H*  l  m  *  **  V  * 


+  -k  C  o 

E  (Rj)  «  N  ^  L  e  1 


&  ,  • 

nk 


Since  this  transformation  is  unitary,  the  ground  state  of  a  one-electron  system 
could  be  written  as  a  determinant  of  Wannier  functions.  This  state  could  include 
the  self-consistent  band  structure  in  a  one-electron  formulation.  The  polarizability 


(6) 


between  the  non-interacting  electron-hole  pair  is  then 


N12,34  (u>) 


|  <1.4  *2,3 


(f  -  f_) 


n 


E1  -  E2  •  w  -  i° 


(3) 


1*  n^,i^,  etc*  *  * . 

where  f^is  the  occupation  number  for  band  n.  Thus  we  have  states  1  and  2  as  a 

particle  and  hole  or  vice-versa,  and  the  poles  of  N°  are  just  the  excitation 

energies  between  occupied  and  unoccupied  bands*  Since  the  dielectric  matrix  is 

£  (w)  B  1  -  V  N°  (oj)  (A) 

the  poles  of  N°  are  also  the  poles  of  £.* 

20 

Hanke  and  Sham  have  shown  how  to  include  particle-hole  interactions  in  this 
matrix  notations*  Their  formula  for  the  polarizability  is 

g,  -  N°  (l+^N0)"1  (5) 

where  is  the  matrix  of  exchange  interactions  between  the  particles  and  holes* 
The  sum  of  all  polarizations  is  then 

£  “  £  [i-  X  !i]  ~1  (6) 

and  the  total  inverses  dielectric  function  is 

fc"1  M*1  +  VS  (w).  (7) 

Note  that  if  JJ  «♦  •»,  S_->  -V-^  and  -*•  0. 

Mow  if  long-range  screening  is  sizable ,  as  is  the  case  in  many  metals,  then  V 

will  dominate  since  the  complete  exchange  term  should  be  screened  (i*e.,£  1  vf ) , 

and  the  resulting  poles  of  e  ^  will  be  just  the  plasmon  poles*  On  the  other  hand, 

in  Insulating  systems  and  in  many  semiconductors,  the  term  is  important  and 

screening  is  not  complete*  Then  the  pole  structure  will  be  different* 

To  Investigate  the  pole  structure  when  is  important,  let  us  consider  the 

simple  case  in  which  diagonal  elements  of  vf  and  V  are  dominant*  Thus  we  have 

V*  v 

Vph,ph 

V*-  V 

PPfhh 


(7) 


where  p(h)  refer  to  particle  (hole)  states,  and  p  and  h  can  be  on  different  sites. 
Now  if  we  rewrite  in  Eq.  (6)  as 


-i. 

S(u>)-  [  N%)  -  (V  -  H  Vs)] 


then  since  all  matrices  are  now  diagonal,  the  poles  in  S  are  simply 

w*E-El-V  ,.+2Vt  , 

p  h  PP,hh  ph,ph  (o) 

after  summing  over  all  spin  states.  This  is  just  the  singlet  excitation  energy 

from  state  h  to  state  p  in  first  order.  Thus  in  this  case,  the  poles  £  ^(u>) 

are  just  shifted  from  the  poles  of  N^0,  but  they  still  correspond  f.o  i  tation 

energies.  This  result  is  not  too  surprising  when  one  recalls  that  fr  Aly  £ 

must  give  the  entire  one-particle  excitation  spectrum,  including  pla*  and 

21 

particle  -hole  excitations 

C.  A  Tight  -  Binding  Approach  to  Calculations  of  Excitation  Structure 
The  above  example  shows  that  there  are  cases  where  the  contention  of 
Inkson  and  Anderson  does  not  apply  and  particle  -  hole  excitations  are  important. 
Indeed  there  are  many  systems  besides  insulators  where  this  is  true.  Of  interest 
in  this  work  are  the  transition  metals  and  also  Cu.  As  one  proceeds  across  the 
transition  metal  series,  the  d-bands  became  narrower  and  particle-hole  excitations 
become  quite  important  so  that,  for  example,  Ni  has  a  complicated  electron  energy 
loss  spectrum.  This  is  true  for  transition  metal  compounds  as  well,  and  Cu^O  has, 
of  course,  strong  exclton  contributions  in  the  loss  spectrum.  In  most  of  these 
cases  the  strongest  excitations  are  localised,  and  local  field  effects  and 
symmetry-breaking  due  to  creation  of  the  hole  are  significant.  Thus  one  must  in¬ 
clude  off-diagonal  matrix  elements  in  any  formulation  of  £  This  is  a  completely 
different  situation  than  the  one  considered  by  Inkson  and  Anderson  who  showed  that 
off-diagonal  elements  were  not  important  in  their  case  •  No  one  has  considered 
in  detail  the  contributions  of  the  off-diagonal  elements  of  £  ^  to  the  effective 
interaction. 


(8) 


We  wish  a  model  for  local  field  effects  which  takes  into  account  the  ground 
state  band  structure  of  our  system.  However,  most  calculations  of  ef1  which  include 
band  structure  effects  are  based  on  orbitals  derived  from  the  Bloch  states  of  the 
band  structure.  These  calculations  are  quite  complicated  and  expensive  which  is  one 
reason  why  no  detailed  studies  of  local  field  effects  have  been  made.  Further¬ 
more,  these  calculations  can  only  be  made  for  a  limited  number  of  perfect  crystals, 
and  we  would  want  a  formulation  that  could  study  trends  in  several  systems  so  that 
we  could  find  the  best  candidates  for  an  excitonic  mechanism.  What  we  need  is  a 
parametrized  model  with  which  we  could  derive  most  of  the  components  of  the  ex¬ 
citation  spectrum  so  that  its  results  would  be  believable  and  which  is  simple 
enough  to  apply  to  a  lot  of  different  structures  and  compositions. 

A  simple  approach  to  this  problem  is  the  tight-binding  theory  used  to  study 

the  electronic  structure  of  transition  metal  alloys.  This  theory  has  given 

22 

results  for  the  density  of  states  in  Cu-Ni  alloys  across  the  whole  range  of 

compositions  which  are  consistent  with  photoemission  and  other  experimental  results, 

and  the  method  has  also  been  used  to  give  detailed  information  on  surface  segregation 

•23 

in  transition  metal  systems.  A  similar  approach  has  already  been  used  by  Onodera 
24 

and  Toyozowa  in  their  important  study  of  excltons  in  mixed  ionic  crystals.  The 
theory  should  be  useful  especially  for  systems  with  narrow  d-bands,  and  it  evidently 
has  the  capability  of  handling  complicated  mixtures  of  components. 

The  d-electrons  in  the  tight-binding  approach  are  described  by  the  usual 
Hamiltonian 

"■  l  +  Jj c  *jci+cj  »> 

Where  includes  all  intra-atomic  Coulomb  interactions  on  site  i  and  ty  is 
the  hopping  Integral  describing  transitions  between  states  1  and  j.  In  the  resulting 
densities  of  states  the  terms  determine  the  positions  of  the  bands  and  the  t^ 
terms  determine  the  widths.  One  can  generalize  this  formalism  straightforwardly  to 
multiple  bands. 
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We  have  used  a  continued-fraction  formulation  of  the  tight-binding  theory 


to  calculate  the  Green’s  function  and  electronic  density  of  states*  In  the 
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continued-fraction  method,  the  local  electronic  Green’s  function  is 


Gii(w)  *  (u>  -  eA  -  A^) 


where  the  self-energy  terms  are 


“"ej  ‘  Aj 


+  0(0, 


+  o(t3) 


- 


and  so  on.  Typically,  one  truncates  the  hierarchy  by  requiring  A^«  A^ •  The 

density  of  states  is  then  derived  from 

D  (co)  *  -JL  I*  (12) 

n 

The  result  in  the  spectral  representation  is  then 

Gii(h,)  “  f  Di(u)  dw  *  <13) 

(*>  -  p 

where  the  usual  pole  structure  is  understood. 

Now  we  may  use  the  results  for  to  calculate  an  excitation  spectrum. 

Using  G  from  Eq.  (13)  we  arrive  at  the  result  for  the  non-interacting  polariz¬ 


ability 


Ni1.ii  *  ^(fi  ~  fi)rAi(u)Ai(u)  dudy* 

"  U  -  U  -  tal 


Using  Eqs»(3)  and  (A),  we  can  then  compute  the  matrix  £  with  the 
electron-hole  Coulomb  and  exchange  interaction  Integrals  given  by 

!  '  J  S*  ■  “ij.kl  ■  i  u*J,kl  •  <151 

The  choice  of  parameters  in  Eq*  (15)  determines  the  excitation  structure*  As 
an  example,  we  consider  again  the  diagonal  interaction  (u  -  A  “  >ij,ki  -  -«• 
Then,  as  before,  Che  Imaginary  part  of  the  sum  of  all  polarizations  is 


(10) 
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Fig. la  Topmost  valence  and  lowest  conduction  band  densities  of 


states  for  C^O.  The  widths  are  the  calculated  widths  (Ref.  25), 
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Fig. lb  The  real  and  Imaginary  parts  of  the  free  particle  polar- 
o  o” 

laatlon  M  .  Note  that  N  Is  aero  below  2.17  eV. 
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If  x  is  zero,  the  excitation  spectrum,  given  by  Im(S(u>)),  is  just  the  convolu¬ 
tion  of  electron  and  hole  states.  If  x  is  within  the  appropriate  range  such 
that 

O' 

«  (o>  >  -  -1  , 

x 

N°"(o.x)  -  0  , 

then  an  isolated  pole  in  the  band  gap  results.  This  is  an  exclton. 

For  our  initial  calculations,  we  have  chosen  a  simplified  exclton  inter¬ 
action.  As  in  the  example  in  Sec.  II.  B,  we  choose 

2,-hf-  -  ~  %  uljf  ij)  Vj.Ki  ( 

m  u  -  h  u 

Then  the  matrix  £  given  in  Eq.  (6)  is  diagonal,  and  we  get 


S(«)-  N°(q)) _ 

l-(u  -  h  u*)N°(o.) 

Now  to  study  the  attractive  part  of  the  local  field  effect,  we  must  look  at 


Rc{e  («)}-  1  +  u  Re{S (to)}. 


From  Eq.  (9)  we  get 


Re  {S}-  VT  Ll-(u  -  h  u  W 


J  (u  -  H  uX)(N°  )2, 


where 


o  o* 

N°  -  N  + 

iM0" 

r  r  i 

i 2  r 

D- 

[  l-(u  -  %  uX)No 

J  +  l 

(U 

and  we  have  suppressed 

the  u>  dependence. 

Since 

X 

u 

larger  than  u,  we  have  (u  -  ux)  <0;  and  thus  Re{S}  will  be  less  than  zero  only 

for  H  <.(u  -  u  )  <0.  Therefore  as  uX  becomes  larger,  we  have  a  larger 

frequency  range  over  which  Re{S)  <0.  The  form  of  Eqs.  (19 »  and  (20)  and  the 

o“ 

era  Involved,  dictate  that  N  must  be  small  (preferably 


m  ml  nm*!  i  <"4* 


zero)  for  R{c  <0.  Therefore  true  excitons  in  Che  band  gap  will  dominace  the 
attractive  interaction.  However,  note  that  all  these  restrictions  may  not  apply  to 
the  case  when  (V  -  %  V*)  is  a  matrix. 

Our  initial  calculations  were  for  Cu^O.  The  lowest  conduction  band  and  highest 
•/«)»#•<#*  I »i mi  are  bull'  foi'Meil  ffmii  the  Cu  Jd  ufbitaJs  accoiditlg  to  band  nltuctufe 

calculations^.  The  exciton  spectrum  is  quite  rich,  but  we  choose  here  only  the 
initial  (n  -  1)  line  of  the  yellow  series27.  The  calculated  densities  of  states 
for  the  valence  and  conduction  bands  are  shown  in  Fig.  1.  These  have  the  widths 
and  relative  positions  of  the  actual  bands.  In  Fig.  2  we  depict  the  results  for 
for  a  range  of  values  of  u  -  u*  appropriate  to  the  exciton  series.  As 
might  be  expected.  Re  {e  }  becomes  more  negative  as  we  increase  the  size  of  u  and  u  . 

We  are  now  extending  this  method  to  realistic  valence  and  conduction  band 
densities  of  states.  We  are  also  including  the  rest  of  the  d-electron  interaction 
integrals  so  that  we  have  a  full  matrix  formulation.  The  result  should  be  a 
believable  picture  of  the  exciton  contribution  to  the  dielectric  response  of  the 
system.  We  are  also  going  to  repeat  these  calculations  with  the  parameters  appropriat 
for  CuCl  and  CdS  systems. 
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Ill,  Development  of  Calculation  Methods 
A*  Molecular  Cluster  Methods 

We  have  Implemented  two  methods  for  calculating  the  self-consistent  electronic 
structure  of  molecular  clusters,  ab-initio  ttartree  Fbck  and  linear-comb ination-of- 
atomic-orbitals-Xa(LCAOXa) •  These  methods  give  a  detailed  description  of  local 
electronic  structure,  and  the  self-consistency  allows  a  description  of  electronic 
rearrangement  even  for  changes  in  complicated  geometries.  In  the  present  case  the 
methods  are  useful  for  studying  the  effects  of  localized  defects  or  bound  electron- 
hole  pairs  on  the  optical  absorption.  The  cluster  approach  can  also  be  used  to 
study  the  contribution  of  localized  excitations  to  the  dynamical  response  of  our 
bulk  systems. 

Although  we  have  done  several  calculations  with  the  Hartree-Fock  method,  we 

will  use  our  Xa  codes  more  extensively  in  the  future.  One  reason  is  that  the  UHF 

code  is  much  more  expensive  than  the  Xa  code  for  a  given  system  owing  to  the 

difference  between  the  detailed  non-local  exchange  calculation  and  the  local  density 

formulation.  Another  reason  is  that  the  Xa  method  usually  gives  accurate  results 

when  compared  with  experiment,  at  least  for  ground  state  properties.  The  Xa 
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method  we  use  is  due  to  Dunlap  et  al,  and  has  the  advantage  over  multiple- 
scattering  formulations  that  it  does  not  use  the  muffln-tim  approximation.  The 
molecular  orbitals  are  in  the  form  of  linear  combinations  of  Gaussian  orbitals 
making  easier  the  calculation  of  matrix  elements  for  system  properties.  The 
code  can  compute  integrals  for  angular  momentum  up  through  f  and  g  so  that  one 
can  study  4d  transition  metals  and  even  rare  earths.  However,  one  would  have  to 
add  relativistic  effects  in  these  cases. 

We  have  modified  the  codes  to  better  approximate  the  bulk  environment.  To 
do  this  we  employ  two  methods.  One  is  to  remove  unpaired  valence  electrons  and  re¬ 
place  them  with  a  compensating  sphere  of  charge  in  which  the  cluster  is  inscribed. 
The  other  method  is  to  include  the  Interactions  with  the  ions  outside  the  cluster 
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by  adding  Madelung  terms*  Both  methods  should  be  useful  for  semiconductor 
systems  like  Cu^O. 

B.  Band  Structure  Methods 

We  have  Implemented  Xa  codes  which  use  an  LCAO  basis  and  thus  are  com- 
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patible  with  our  cluster  methods.  The  codes  were  developed  by  Wang  and  Callaway 
and  include  d-  and  f-  functions  so  they  can  be  used  on  transition  metal  systems. 
However,  the  programs  were  originally  written  only  for  simple,  body-centered  and 
face-centered  cubic  structures  with  one  atom  per  unit  cell.  We  have  now  extended 
them  to  rocksalt  structures  so  that  we  can  look  at  transition  metal  oxides.  In 
the  future  we  will  extend  them  to  more  general  cubic  structures  with  two  atoms 
per  unit  cell  so  that,  for  example,  we  can  look  at  CuCl.  The  band  structure  codes 
are  also  spin-unrestricted,  as  are  the  cluster  codes,  and  we  can  therefore  study 
spin  polarization  effects  which  are  important  in  some  transition  metal  systems. 

With  the  results  of  the  band  structure  calculations,  we  can  look  at  such 
things  as  optical  properties  and  use  the  band  vavefunctlons  and  matrix  elements 
as  a  basis  for  studies  of  the  excitation  structure.  We  may  also  use  the  results 
as  a  starting  point  to  obtain  better  parameters  for  our  tight-binding  calculations. 
And  eventually  we  will  use  the  results  to  provide  a  detailed  environmental  po¬ 
tential  in  a  more  complete  bulk  embedding  model  for  our  cluster  calculations. 
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IV,  Electron  Energy  Loss  Structure  of  Clean  and  Oxidized  Titanium 


In  this  section  we  report  the  completion  of  our  work  on  analyzing  con¬ 
tributions  of  lnterband  transitions  and  plasmons  to  the  excitation  spectra  of 
transition  metals. 

There  is  a  qualitative  disagreement  among  the  experimental  loss  spectra* 
and  in  addition  the  results  have  been  given  different  interpretations  so  that 
very  different  energies  for  plasmons  and  interband  transitions  have  been  re¬ 
ported,  The  discrepancies  in  measurements  on  transition  metalr.  make  it  unclear 
whether  the  plasmon  energy  position  is  given  by  free-electron  theory  or  whether, 
for  example* effects  due  to  the  periodic  lattice  potential  are  important. 

The  characteristic  loss  structure  measurements  fall  into  two  classes: 
measurements  of  losses  suffered  by  electrons  with  energies  of  many  kilovolts  as 
they  are  transmitted  through  thin  films  and  measurements  of  losses  suffered  by 
electrons  with  energies  around  lOOeV  as  they  are  reflected  from  the  material. 

The  former  experiments  are  dominated  by  electrons  which  do  not  undergo  elastic 
scattering  while  the  latter  are  more  highly  resolved  but  contain  electrons  which 
have  been  elastically  scattered  at  least  once.  A  typical  reflection  spectrum  for 

clean  Ti  is  shown  in  the  lowest  curve  (E  -  200eV)  of  Fig,  3.  The  transmission 

P 
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spectrum  differs  in  that  the  main  peak  is  shifted  several  volts  higher  in  energy. 

A  similar  situation  occurs  when  we  compare  measurements  we  have  made  on  Fe/FeO 
systems  with  transmission  experiments.  To  clear  up  the  confusion  and  also  to 
provide  a  basis  for  use  of  EELS  as  a  surface  analysis  tool*  we  made  several 
measurements  of  the  loss  spectra  for  Ti  systems. 

The  analysis  of  the  clean  and  oxidized  Ti  data  in  Figs,  3  and  4  was  given 
in  our  last  report.  In  sunmary*  for  the  clean  TI  spectra  in  Fig.  3*  the  peak 
around  5eV  is  a  surface  plasmon  which  very  quickly  diminishes  with  oxygen  coverage 
or  Increasing  primary  energy.  The  broad  main  peak  Is  the  bulk  plasmon  with 
secondary  surface  and  bulk  excitations  in  its  tall.  The  peaks  at  about  42eV* 


(15) 


60 


50  40  30 

ENERGY  LOSS  ( 


20  10  0 
EV  ) 


Mf»3  A  comparison  of  dean  Ti  loan  spectra  for  different 
primary  energies. 


which  increases  in  intensity  relative  to  the  main  peak  as  primary  energy  in¬ 
creases,  is  mysterious.  It  is  in  the  wrong  position  for  an  interband  excitation 
and  does  not  behave  like  a  collective  mode.  It  may  be  due  to  a  collective 
excitation  triggered  by  an  excited  Ti3p  electron.  For  the  oxidized  Ti  spectra 
in  Fig.  4,  the  lower  curve  shows  the  surface  peak  greatly  reduced  and  the 
plasmon  peak  shifted  to  somewhat  lower  energy.  The  02s  peak  is  now  prominent  at 
about  20eV,  and  the  broad  42eV  peak  persists.  The  upper  curve  for  EP*1000eV 
shows  a  significantly  different  spectrum  with  a  broadened  plasmon  peak  now  shifted 

to  overlap  the  02s  peak.  If  the  shifting  were  to  continue  for  higher  primary 
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energies,  the  spectrum  would  look  much  like  the  transmission  spectrum.  The 
shifting  of  the  bulk  plasmon  peak  with  primary  energy  is  seen  both  in  the  clean 
and  oxidized  systems  and  is  also  seen  for  other  materials. 

What  could  cause  the  movement  of  the  plasmon  peak  to  higher  energies  as 
the  incident  energy  is  increased  in  the  reflection  experiments  and  also  could 
cause  the  difference  between  the  positions  of  the  plasmon  peak  in  the  reflection 
and  transmission  experiments?  A  possible  explanation  of  this  movement  is  that  it 
is  caused  by  plasmon  dispersion.  The  RPA  formula  for  dispersion  is 


so  that  non-zero  momentum  transfer  would  always  result  in  a  positive  energy  shift. 
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Now  in  the  transmission  experiments,  the  acceptance  angle  of  the  detector  is 
5  nr  ad.  so  that  for  Incident  energies  of  35  keV  the  momentum  spread  for  the 
electrons  transmitted  in  the  forward  direction  is  well  below  the  Fermi  momentum. 
However,  if  there  is  enough  scattering  of  the  electrons  within  the  sample,  a 
significant  number  of  the  electrons  undergoing  plasmon  losses  and  then  entering 
the  detector  would  transfer  considerable  momentum  to  the  platoons.  We  have 
applied  to  the  electrons  a  model  derived  for  multiple  scattering  of  ions  in  a 
screened  Coulomb  region^  and  have  found  that  the  electron  beam  does  spread  enough 
to  allow  a  significant  number  of  large  momentum  plasmon  losses  to  be  detected. 
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Pig.4  Lot*  spectra  of -final  oxidized  surface  for  EP*200  eV  and  EP-1000  eV 
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This  effect  would  cause  a  shift  to  higher  energy  and  a  broadening  of  the  plasmon 

peak.  A  similar  effect  would  also  take  place  in  the  reflection  experiments 

since  as  the  incident  energy  increased  the  electrons  would  penetrate  deeper  into 

the  sample  and  be  scattered  more.  This  would  explain  the  dependence  of  the 

plasmon  peak  position  on  incident  energy. 

There  are  evidently  no  measurements  of  the  plasmon  dispersion  for  Tlt  but 
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measurements  for  AJt  indicate  a  dispersion  of  as  much  as  13  eV.  If  a  plasmon 
dispersion  of  several  eV  with  a  moderate  increase  in  peak  width  is  possible  for 
Ti,  this  effect  could  explain  the  discrepancy  between  the  transmission  and  re¬ 
flection  measurements.  In  that  case  the  reflection  peak  would  be  closer  to  the 
true  zero-momentum  volume  plasmon  energy.  Before  any  definite  conclusion  can  be 
reached,  however,  one  must  do  a  calculation  with  a  more  realistic  model  of  electron 
scattering. 

These  results  resolve  some  of  the  conflicts  between  transmission  and  re¬ 
flection  loss  spectra.  They  also  reveal  that  care  must  be  exercised  in  taking 
and  analyzing  the  data.  In  particular  for  a  reactive  surface  such  as  that  of 
Ti,  oxygen  adsorption  can  greatly  change  the  appearance  of  the  spectrum.  Our 
analysis  provides  a  basis  for  understanding  the  origins  of  the  peaks  in  the 
spectra,  how  they  change  as  conditions  are  varied,  and  especially  what  happens 
to  the  surface  peak. 
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1.  "Optical  Constants  of  TiO^".  to  be  published  in  the  Handbook  of 
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2.  "Local  Fields  due  to  Excitons  and  Their  Effect  on  Superconductivity". 
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Mr.  David  Luedtke  -  Ph.D  candidate. 
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VIII «  Interactions 


We  had  the  following  interactions  in  connection  with  this  research  during 

the  period  October  1,  1980  to  September  30,  1981: 

1*  With  Dr*  C.  S*  Wang  at  the  University  of  Maryland  on  the  modification 
a 

•  and  use  of  LCAO-Xa  band  structure  codes  developed  by  her. 

2.  With  Dr.  Brett  Dunlap  of  N.R.L.  about  using  and  modifying  LCAO-Xa 
cluster  codes  he  has  developed. 
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